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We show that if the sources of ultra-high energy cosmic rays (UHECRs) with energies E > 10 19 eV 
are associated with galaxies, then the production of UHECRs must occur in transient events. Our 
galaxy is currently at a dim state in between transients. The time interval between transients is 
> 10 4 ' 5 min{l, (7min/10 3 ) -0 ' 6 } yr for 7 m m ^ 10 7 , where 7 m in is the minimum Lorentz factor to 
which protons are accelerated in the transients. This constraint is satisfied by 7-ray bursts. 



PACS numbers: 98.70.Sa, 96.40.-z, 98.62.-g 

The origin of the observed cosmic rays at different en- 
ergies is still unknown (see for recent reviews). As 
illustrated in Fig. 1, the cosmic ray spectrum changes its 
qualitative behavior as a function of particle energy; it 
steepens around ~ 5 x 10 15 eV (the "knee") and flattens 
around 5 x 10 18 eV (the "ankle"). Below - 10 15 eV, the 
cosmic rays are thought to originate from Galactic su- 
pernovae. The composition is dominated by protons at 
the lowest energies, and the fraction of heavy nuclei in- 
creases with energy. The proton fraction at ~ 10 15 eV is 
reduced to ~ 15% [|]||]. At yet higher energies, there 
is evidence that the fraction of light nuclei increases, 
and that the cosmic-ray flux above 5 x 10 18 eV is again 
dominated by protons |(|. This composition change and 
the flattening of the spectrum around 10 19 eV (see Fig. 
2) suggest that the flux above and below this energy 
is dominated by different sources. Since the Galactic 
magnetic field can not confine protons above 10 19 eV, 
it is believed that the nearly isotropic cosmic ray flux at 
E > 5 x 10 18 eV, originates from extragalactic sources. 
The small, but statistically significant, enhancement of 
the flux at E < 3 x 10 18 eV near the Galactic plane 
suggests a Galactic origin at these lower energies. 

The Milky- Way disk shows prominently relative to 
the extragalactic background in electromagnetic radia- 
tion ranging from radio to X-ray wavelengths. This is 
a consequence of the fact that our location within the 
Galaxy is not representative of a random point in the 
universe, where the observed radiation would be nearly 
isotropic. The Galactic prominence must also apply to 
UHECRs with energies > 10 19 eV. At these energies the 
Galactic magnetic field is unable to isotropize the par- 
ticle orbits, and the Galactic disk should be on aver- 
age brighter than the extragalactic background by orders 
of magnitude. Contrary to this expectation, the actual 
anisotropy of the observed UHECR is < 4% in the di- 
rection of the Galactic disk |R|P]. In this paper, we 
therefore argue that any generic galactic production of 
UHECRs must occur in transient events with a short 
duty cycle. Using existing data for the spectrum and the 
confinement time of cosmic rays at lower energies, we set 
limits on the event rate and the minimum Lorentz factor 



of the accelerated particles in these transients. 
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FIG. 1. Data points from [2], showing the differential flux 
of all cosmic-ray particles multiplied by E 3 . The dotted line 
is a power-law fit to the observed proton flux [3]. The solid 
line shows the expected volume- averaged, extragalactic cos- 
mic-ray flux based on the observed UHECR flux (see text 
and Fig. 2). The corresponding minimum, time-averaged, 
Galactic flux is shown as the dashed line. 

The extragalactic cosmic-ray flux. In Fig. 2 we com- 
pare the UHECR spectrum measured by the three exper- 
iments with the largest current exposure, Yakutsk, Fly's 
Eye and AGASA, with that expected from a homoge- 
neous cosmological distribution of sources, each gener- 
ating a power-law, differential number flux of high en- 
ergy protons, dn/dE oc E~ 2 . This power-law slope is 
expected for astrophysical sources which accelerate par- 
ticles in strong collisionless shocks [OlO], and is found 
in supernovae fHI and 7-ray bursts |fl6| . Experimental 
differences in the absolute flux calibration at 3 x 10 18 eV 
yield a systematic over-estimation of the event energies 
by ~ 20% (~ 10%) in the Yakutsk (AGASA) experiment 
as compared to the Fly's Eye experiment; in Fig. 2 we use 
the Yakutsk energy normalization. The calculation of the 
model flux follows ref. fTif] , and assumes a flat universe 
with fi m = 0.3, fl A = 0.7 and H = 65 km s" 1 Mpc" 1 . 
The generation rate of cosmic rays per unit comoving 
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FIG. 2. The UHECR flux expected in a cosmological 
model for which high-energy protons are produced at the rate 
given by Eq. (1) (solid line) compared to the Fly's Eye [9], 
Yakutsk [10] and AGASA [11] data (adopted from [12]). The 
flux error bars correspond to la. The highest energy points 
are derived assuming that the detected events represent a uni- 
form flux over the energy range 10 20 eV - 3 x 10 20 eV. The 
dashed line is the sum of the cosmological model flux and 
the Fly's Eye fit to the Galactic heavy nuclei component, 
J oc E~ 3 - 5 [9]. 



volume is assumed to trace the redshift evolution of the 
luminosity density of bright quasars [Q, which is also 
similar to the evolution of the cosmic star formation rate 
iP^t , namely ucr(z) oc (1 + z) a with a « 3 at low red- 
shifts z < 1.9, ticr(z) — const for intermediate redshifts 
1.9 < z < 2.7, and an exponential decay at high redshifts 
z > 2.7. The cosmic-ray spectrum at energies > 10 19 eV 
is only weakly dependent on the cosmological parameters 
or the assumed redshift evolution, because the particles 
at these energies originate from distances shorter than a 
few hundred Mpc. The spectrum and flux at E > 10 19 eV 
is mainly determined by the present [z = 0) generation 
rate and spectrum, which for the model shown in Fig. 2 
is given by 
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The suppression of the flux above 10 197 eV is caused 
by the energy loss of high energy protons in their inter- 
action with the microwave background, i.e. the "GZK 
cutoff" (^0|. The available data does not allow to de- 
termine the existence (or absence) of the "cutoff" with 
high confidence. The AGASA results show an excess 
of events (at a ~ 2.5a confidence level) compared to 
the model prediction above 10 20 eV. This excess is not 
confirmed, however, by other experiments; preliminary 
results from the new HiRes experiment are consistent 
with the Fly's Eye data which does not show this excess 
pjj . Since the flux at 10 20 eV is dominated by sources 
at distances < 100 Mpc over which the distribution of 



known astrophysical systems (such as galaxies or galaxy 
clusters) is inhomogeneous, significant deviations are ex- 
pected at this energy relative to the model predictions 
presented in Fig. 2 under the assumption of a uniform 
source distribution p7 22|. 

Galactic enhancement of observed flux. If the produc- 
tion rate of high energy protons given in Eq. (|l|) extends 
to lower energies, then a universal flux of extragalactic 
protons would be produced as shown by the solid line in 
Fig. 1. If extragalactic cosmic rays originate in galaxies 
similar to our own, then the average Galactic flux would 
be enhanced relative to this extragalactic background by 
two factors. The first factor is cosmological and is inde- 
pendent of the confinement of cosmic rays by the Galactic 
magnetic field. Ignoring confinement for the moment, we 
expect that the extragalactic flux would be diluted rel- 
ative to the Galactic flux by the volume filling factor of 
galaxies tigVg (where nc is the number density of galax- 
ies and Vq is the volume of a galaxy), and be amplified by 
the ratio between the emission time (which is of order the 
Hubble time, fn) and the escape time from the Galactic 
disk (which is of order the vertical crossing time of the 
scale-height of the Galactic disk, h/c). The combination 
of these ratios yields a Galactic enhancement factor of 
the cosmic ray flux J(E), 
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where Aq = Vq/H is the typical cross-sectional area of 
a galactic disk, and the subscripts G and C are used to 
denote the Galactic and cosmological values. The numer- 
ical value of /i ~ 10 3 is obtained for the characteristic 
Aq ~ 7r(5 kpc) 2 and uq ~ 3 x 10 -3 Mpc" 3 of L+ galaxies 
similar to the Milky- Way [§3| . 

A second enhancement factor originates from the con- 
finement of cosmic rays by the Galactic magnetic field. 
The time that a cosmic ray spends in the disk is larger 
than h/c by a factor /a- This factor is directly deter- 
mined by observations of daughter nuclei produced in 
interactions of cosmic ray nuclei with H and He nuclei 
in the Galactic disk. The abundance ratio of daugh- 
ter to parent nuclei determines the average grammage 
E (defined as the path length integral of the gas den- 
sity) that cosmic ray nuclei traverse prior to their es- 
cape. The factor fa is given by the ratio S/Sdisk, where 
Edisk = 2 x 10~ 3 g cm~ 3 is the surface mass density of 
the Galactic disk. For cosmic rays with energies of up to 
1 TeV per nucleon jgl 
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where Z is the atomic charge and a ~ 0.6. The ex- 
trapolation of this scaling relation to a proton energy 
E r~j 10 16 eV gives S ~ Edisk or w 1, which implies 
that the extrapolation can not be extended to yet higher 
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proton energies since S must be larger or equal to Sdisk- 
Note that the extragalactic background is also enhanced 
by a factor f% if measured within the Milky- Way galaxy. 
Hence the ratio between the local values of the mean 
Galactic and extragalactic fluxes remains /i even in the 
presence of confinement. 

Figure f shows that the observed energy flux E 2 J{E) 
at 10 15 eV is larger by a factor of ~ 10 3 than its value at 
10 19 eV. This factor is coincidentally comparable to the 
enhancement factor f\. If, as commonly assumed, the 
cosmic rays at 10 15 eV are Galactic in origin while those 
at 10 19 eV are extragalactic, then the energy production 
rate per galaxy at these two cosmic-ray energies must be 
comparable. Since the particle spectrum expected for ac- 
celeration in astrophysical shocks, E 2 dN/dE = const, 
implies equal amount of energy per logarithmic energy 
interval, the above coincidence may seem to suggest that 
both the 10 15 eV and 10 19 eV cosmic rays are produced 
by the same source population [p5f . However, the asso- 
ciated Galactic component of cosmic rays with energies 
between 10 15 eV and 10 19 eV appears to be missing. The 
spectrum steepens instead of maintaining a J(E) oc E~ 2 
slope in Fig. 1, as expected if the same sources were 
producing the observed cosmic-ray flux at 10 15 eV and 
10 19 eV. The lack of strong anisotropy in the direction 
of the Galactic disk at 10 19 eV supports this inference. 
Where are the Galactic counterparts to the extragalac- 
tic 10 19 eV cosmic rays? To reconcile the data in Fig. 
1 with the enhancement factors in Eqs. (||) and (|^), we 
argue that the UHECRs must be produced in transient 
events and that our galaxy is currently at a dim state in 
between transients. 

Constraints on transient properties. The minimum 
flux of lower energy cosmic rays produced by the UHECR 
source population corresponds to an extrapolation with 
a spectral slope of J(E) oc E~ 2 , which is the shallowest 
slope possible in shock acceleration (analysis of the ob- 
served spectrum at E > 10 19 eV implies that the source 
spectral slope can not be shallower than J(E) oc 
fin\). The dashed line in Fig. 1 shows this extrapolation 
together with the expected enhancement factors /i and 
/2 from Eqs. (||) and (||), assuming conservatively that 
Eq. (||) with a = 0.6 holds for protons up to E ~ 10 16 eV. 
At energies exceeding ~ 1 TeV, the expected mean Galac- 
tic flux of protons exceeds the observed value. The fact 
that the bulk of the Galactic cosmic rays in this energy 
range are missing is supported by the isotropy of the ob- 
served cosmic rays at 10 19 eV. The Galactic field is un- 
able to isotropize the distribution to the observed level 
since the Larmor radius of a particle of atomic charge Z 
in the Galactic field B is - 10 kpc(£/10 19 eV)/(B/fiG), 
significantly larger than the Galactic scale height of this 
magnetic field. 

Our assumption that Eq. (^|) holds for protons up 
to E ~ I0 16 eV is conservative since it implies that 
£ = Sdisk at E ~ 10 16 eV. If £ is higher, e.g. due to a 



smaller value of a beyond 1 TeV per nucleon, then the 
expected Galactic counterpart to the extragalactic flux 
will be higher than shown in Fig. 1. 

The confinement time of cosmic rays with an energy 
of ~ 1 GeV per nucleon is measured to be ~ 10 75 yr, 
based on the survival fraction of radioactive cosmic-ray 
nuclei ^6|. As the cosmic-ray energy increases, the de- 
cline in £ described in Eq. (Q) may be caused either by 
a decrease in the confinement time or by a decrease in 
the time-averaged gas density through which the cosmic 
rays propagate (e.g., due to an increase in the thickness 
hen of the disk in which the cosmic rays are confined to a 
value larger than the thickness h of the Galactic gaseous 
disk). Assuming, conservatively, that the decrease in £ 
is only due to a decrease in the confinement time, we get 
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In order that the Galactic counterpart to the extra- 
galactic ~ 10 19 eV cosmic ray protons will not exceed 
the observed proton flux above ~ 1 TeV, the time be- 
tween Galactic transients must satisfy Ttrans > /vis x 
min{T conf (E = 7mi n GeV), r conf (E = 1 TeV)}. Here, 7 min 
is the minimum Lorentz factor of the accelerated protons 
and /vi s is the fraction of all Galactic transients which are 
visible to us. Since diffusion perpendicular to the Galac- 
tic disk drains the cosmic-ray flux from a distant source, 
a transient would be visible to us only if it has occurred 
within a distance comparable to the scale height hcR of 
the cosmic-ray disk. The estimated scale-height, based 
on the time-averaged density through which the cosmic 
rays propagate, is ~ 3 kpc Q implying / vis ~ 0.1. We 
therefore obtain 
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as long as 7 m in ^ 10 7 . At higher values of 7 m i n for which 
£ ~ Sdisk and the cosmic rays escape freely from the 
Galactic disk, the only constraints are that Tt ra ns must 
be longer than the exposure time of modern experiments 
(several decades) and that the transients have a short 
duty cycle. 

Radioactive dating provides only a conservative lower 
limit on the confinement time due to the potential exis- 
tence of an extended cosmic-ray halo. We note that the 
above conclusions remain unchanged even if the 10 19 eV 
cosmic rays are heavy nuclei since the factor /i is inde- 
pendent of composition. 

Discussion. Among the most likely astrophysical 
sources for the unsteady production of UHECRs are 
Gamma-Ray Bursts (GRBs) and active galactic nuclei 
(AGNs). While AGNs (such as SgrA* in the Milky Way 
galaxy) are known to be transient just as GRBs, the tran- 
sient properties of GRBs are better-determined and allow 
a more detailed comparison with the required properties 
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of UHECR sources. The GRB rate per comoving volume 
today is ~ 0.5 x 1CT 9 Mpc" 3 yr _1 0. This rate is de- 
rived based on the observed value at a redshift z ~ 1, 
assuming that the rate follows the redshift evolution of 
the cosmic star formation rate and that the emission of 7- 
rays by the bursts is isotropic. If the emission is confined 
to a double sided jet with a small angular radius 9, then 
the burst rate is increased by a factor 4tt/2tt9 2 . GRB ob- 
servations indicate a characteristic opening angle 9 > 0.1 
p8| , p9[ . The implied temporal separation between bursts 
in L+ galaxies is thus tgrb ~ 10 4 5 (^/0.1) 2 yr, consistent 
with the constraint in Eq. (|5|). 

Under the assumption that supernovae dominate the 
cosmic-ray production below 10 15 eV, the value of 7 m i n 
for GRBs with tqrb ~ 10 4 ' 5 yr cannot be substantially 
lower than ~ 10 3 . This constraint is naturally satisfied in 
GRBs since the plasma within which particles are accel- 
erated, is inferred to expand with a bulk Lorentz factor 
r ~ 10 2 - 5 (e.g. go)). 

The required energy production rate at 10 19 eV is ~ 
10 44 erg Mpc~ 3 yr -1 implying that the cosmic- 

ray energy per burst is 10 51 (rb U rst/10 4 ' 5 yr) ergs. This 
energy release is surprisingly close to the energy carried 
by 7-rays in GRBs for T tralls = 10 4 ' 5 yr (corresponding 
to 9 r~j 0.1 = 6°) p9[| , implying that the local energy gen- 
eration rate of 7-rays in GRBs is similar to the UHECR 
energy generation rate |3^]Q. This similarity, the agree- 
ment with the transient rate constraint, and the fact that 
protons may be accelerated to energies > 10 20 eV by col- 
lisionless internal shocks of GRBs |33| support the sug- 
gested association between GRBs and UHECR sources 
©0. 

The shocks produced due to converging flows during 
the formation of large scale structure in the intergalactic 
medium are also expected to accelerate cosmic rays [ j35| . 
For the characteristic magnetic field strength and scale 
of these shocks, we estimate that the maximum proton 
energy is ~ 10 17 eV and that the intergalactic cosmic rays 
may approach the observed flux around this energy but 
fall significantly below the observed flux at lower energies. 

In closing, we emphasize that the Earth was exposed 
to a time-averaged UHECR flux that is greater by three 
orders of magnitude than the currently measured value. 
During transients the enhancement is much larger^]. Any 



*The AG AS A excess above 10 20 eV has been used to argue 
that GRBs can not produce the observed cosmic rays above 
10 20 eV j3lj. As explained above, the excess is not seen in 
the Fly's Eye, the Yakutsk and the HiRes data, and if real is 
likely associated with source inhomogeneities. 

^The transient event rate and intensity might depend on 
particle energy. For beamed GRBs, the bursts at > 10 19 eV 
will be less frequent but more intense than those at much 
lower particle energies, since the collimation of UHECRs at 



fossil record of this intense cosmic-ray bombardment (e.g. 
through the detection of stable or radioactive daughter 
nuclei in ancient rocks) would provide an important test 
of our conclusions. 
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